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Floppy membranes are tensionless surfaces without extrin- 
sic stiffness whose fluctuations are governed by fourth-order 
bending elasticity. This suppresses spiky superstructures and 
ensures that floppy membranes remain smooth over any dis- 
tance, with Hausdorff dimension Dh = 2, in contrast to sur- 
faces with stiffness which are rough on the scale of some finite 
persistence length. 
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1. Under deformations, fluid membranes behave ap- 
proximately like ideal tensionless surfaces with curvar- 
ture stiffness , and have a model energy 



H = n J d 2 £ yfg Tr C 2 



(1) 



where C a b is the second fundamental form of the surface 
described by an embedding function x = x k 1 , £°2) . The 
symbol g denotes the determinant of the induced metric 
g a b = <3 a x<9 b x. 

Thermal fluctuations are known to soften the curvature 
stiffness with increasing membrane size, such that there 
exists a finite persistence length C = Co exp(47r«;/3T), 
where £0 is the molecular size, beyond which the mem- 
brane looses its stiffness completely and begins exhibiting 
a surface tension proportional to (~ 2 , which was initially 
absent. Thus the tangential correlation functions have 
only the range £, and surfaces much larger than £ appear 
rough. The typical behavior of bilipid vesicles can there- 
fore be observed in the laboratory only for sizes smaller 
than C- Surfaces much larger than £ will crumple and fill 
the embedding space with spiky structures (JJ^). 

Apart from these common membranes, nature may 
provide us also with another type, not described by the 
Hamiltonian (Q). ff the molecules in a bilayer are strongly 
conical, the bending stiffness can be zero or negative |6| . 
Such a situation can also arise for charged or dipolar 
molecules . In this case we shall speak of floppy mem- 
branes. The purpose of this note is to study the statisti- 
cal properties of such floppy membranes, which will turn 
out to be quite different from those of ordinary mem- 
branes. In particular we shall find that, in contrast to 
ordinary membranes, floppy membranes without tension 
and stiffness are smooth over long distances. If the stiff- 
ness is negative, floppy membranes are able to form disor- 
dered superstructures similar to those recently observed 
in the laboratory ||. These are thought to be molten 
versions of the egg carton-like crystalline arrangement of 



local maxima and minima on the surface. This has not 
yet been confirmed experimentally, but was suggested by 
recent numerical simulations j9|. 

2. In order to describe floppy membranes, we must 
stabilize their fluctuations by adding to the energy (Q) 
a higher-gradient term. Here we shall consider only one 
of several possibilities, focusing our attention upon the 
following Hamiltonian: 



(2) 



H=r J d 2 ^g+- J d 2 ZJgTrC 2 

+ i / [g cd v a c ac v b c bd + TrC 4 ] 



The properties of a surface with this Hamiltonian will be 
studied non-perturbatively in the limit of a large number 
D of embedding dimensions. 

The first term in (||) parametrizes the surface tension; 
the third term provides the surface with the stabiliz- 
ing higher-order bending stiffness, whose parameter m 
has the dimension (energy x surface) -1 . Although this 
term is irrelevant in a perturbative renormalization group 
analysis, it becomes relevant non-perturbatively in the 
limit ID — ► 00 by a mechanism familiar from the three- 
dimensional Gross-Neveu model jl0| . The new term sta- 
bilizes the surface against growing spikes and makes it 
smooth over long distances. 

The Hamiltonian (^) can be reformulated alternatively 
in terms of the tangent vectors 9 a x of the surface, or in 
terms of the normal vectors n: 



+ \ J d 2 ^g n 



-kTT 



(3) 



where K is the Gaussian curvature. This form exposes 
an important physical aspect of the model. As pointed 
out in pT| ] , the first term in the second line is analogous 
to the the continuous version of the Heisenberg model 
of ferromagnets p2| , albeit with an additional integra- 
bility condition for the n(£)-field. It tries to make an 
ordinary membrane with a positive stiffness smooth cor- 
responding to a ferromagnetic alignemcnt of the normal 
vectors. The fact that ordinary membranes cannot be 
smooth over long distances has its parallel in the absence 
of an ordered phase in the two-dimensional Heisenberg 
ferromagnet. By the same analogy, we see that our new 
term introduces (apart from an intrinsic term K 2 ) an 
antifcrromagnetic next-to-nearest neighbours interaction 



1 



between normal vectors. This generates frustration, and 
it is due to this non-local interaction that the surface 
can have an ordered phase after all, although with an 
antifcrromagnetic type of order. 

To exhibit these features analytically, a formulation of 
the Hamiltonian (J2J) in terms of the tangent vectors will 
be most convenient: 



H = \ J d 2 £^g g ab V a x^ (r - nV 2 + ±-V^j V b 2 



(4) 



3. We analyze the model (|J) in the large- D approx- 
imation along the lines of Refs. |i|Jl^|. To this end we 
introduce a Lagrange multiplier matrix A ab to enforce 
the constraint g a b = <9 a xd;,x, extending the Hamiltonian 
(1 to 



#ext - H + ~ J d 2 i^ X ab (9 a x9 h x - g ab ) 



(5) 



Then we parametrize the surface in a Gauss map by 
x(0 - (i = 3,...,D), where -R x /2 < 

Ci < Ri/2, -R2/2 < 6 < R2/2 and <j> l (£) describe the 
(D-2) transverse fluctuations. In the limit of large D, 
to be studied here, the large number of components sup- 
presses the fluctuations of A afc and g a b. These fields take 
extremal values which, for large surface areas, are homo- 
geneous and istropic: g a b — P Sab, A ab = A g ab . Thus we 
may replace @ for large D by 

H cxt = J d 2 ^ [r + X(l-p)} 

+ \ J d 2 H K (-V 2 ) d a ^ , 

where K represents the differential operator 

K (-V 2 ) = r + A - kV 2 + — X> 4 . 

m 



(6) 



(7) 



Integrating out the transverse fluctuations, always for 
large areas, we get the free energy 



F=A cxt [r + A(l-p)] 
D-2 



+A C 



8tt 2 



p d 2 p\n[p 2 K(p 2 )] , 



(8) 



where, for simplicity, we have chosen natural units by 
setting (3 = l/k^T = 1 and A ext = R1R2 is the extrinsic, 
projected area in the coordinate plane. The factor (D — 2) 
in the second term ensures that, for large D, the fields 
A and p are extremal and satisfy thus the saddle-point 
("gap") equations 
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= / (r, k, m, A) , p = — 



m, A) 



(9) 



The prime denotes derivatives with respect to A and the 
saddle-point function / is defined by 



f(r,K,m,X) = A - 



D-2 
8tt 



dp p In [p 2 K (p 2 )] . (10) 



Inserting (|9|) in (^) we get F = (r + A) A cxt showing that 
r p h = (r + A) is the physical surface tension. 

We now introduce the following combinations of the 
parameters of the model: 



i? 2 = - v /m(r + A) + — , I 2 = -v / m(r + A) — . 

(11) 

In terms of these, the kernel V can be written as 

mK (p 2 ) = (R 2 + I 2 ) 2 + 2{R 2 - I 2 ) p 2 +p\ (12) 

from where we deduce the stability condition for the ho- 
mogeneous saddle-point as being R 2 > 0, insuring that 
the spectrum of transverse fluctuations is positive for all 
p > 0. Since we shall mostly consider negative or van- 
ishing stiffnesses n we shall also assume that I 2 > 0, 
so that both R and / are real. Note that the spec- 
trum of transverse fluctuations depends drastically on 
the sign of the stiffness: for negative k we have I > R 
and K develops a minimum at p = V I 2 — R 2 . Corre- 
spondingly (but for a slightly higher value of I/R), the 
spectrum E (p 2 ) = p 2 K (p 2 ) develops a roton-like mini- 
mum, which, for R/I < 1, lies at p ~ / (l - 8i? 2 /3/ 2 ), 
as shown schematically in Fig. 1. 
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FIG. 1. The roton-like minimum in the spectrum of trans- 
verse fluctuations for I 3> R. 

Having established the stability conditions we proceed 
to the evaluation of the saddle- function. This contains an 
ultraviolet divergent integral which must be regularized. 
To this end we use the standard dimensional regulariza- 
tion, computing the integral in (2 — e) dimensions. For 
small e, this leads to 



/ (r, k, m, A) = A + — (R 2 - I 2 ) In 
Ait 



R 2 



A 2 

1 nT {* I 2 -R 2 ' 

Rl — h arctan 

2tt I 2 2RI 



(13) 
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where A = p exp(2/e) and p is a reference scale which 
must be introduced for dimensional reasons. The scale A 
plays the role of an ultraviolet cutoff, diverging for e — > 0. 

In order to distinguish the various phases of our model 
we compute two correlation functions. First, we con- 
sider the orientational correlation function g a b{£, — £,') = 
(da^iO db^iC)} f° r the normal components of tangent 
vectors to the surface. From (||) this is given by 



d 2 P 



1 



K(p 2 ) 



(14) 



In terms of R and /, the Fourier components can be 
writtens as 



1 



K(p 2 ) 2RI 



Im 



1 



P 



(R-il) 2 ' 



(15) 



from where we obtain the following exact result for the 
diagonal elements g = g aa of (114) : 



2ng{d) 



2RI 



Im K ((R— il) y/pd) , 



(16) 



where d = \$— £'\ and Kq is a Bcsscl function of imaginary 
argument E4]. 

Secondly, we compute the scaling law of the distance 
6e in embedding space between two points on the surface 
when changing its projection d on the reference plane. 
The exact relation between the two lengths is 



d 2 E = d 2 +j2(W(o-^(a" 



(17) 



exhibiting short-range orientational order. For short 
distances d <C 1/Ry/p, the surfaces behave as two- 
dimensional objects. If p becomes large we have a region 
1/Ry/p <C d <C l/^/27rr p h in which d E scales logarith- 
mically with d and distances along the surface become 
large. The transition to this regime happens on the scale 
of the persistence length d E L — 1/ */f^h. Above this scale 
world-sheets are crumpled, with no orientational corre- 
lations (if the tension is not large enough to dominate 
over the entire surface, causing p ~ 1). This phase corre- 
sponds to the familiar behaviour of stiff membranes [Q . 

For large negative stiffness k, in contrast, we have / 3> 
R, the asymptotic region is d » 1/Iy/p, and R can be 
neglected in Kq for d <C 1/ R\fp- In this region we have 



9(d) 



8RI 



Jo (Iy/pd) 



(19) 



with 1/Ry/p playing the role of an infrared cutoff for 
the oscillations on the scale 1 / lyfp over which the Bessel 
function Jo varies. We have thus a new scale (in embed- 
ding space) 



d° E 



32Rr ph 



(20) 



on which the transverse fluctuations create oscilla- 
tions characterized by the "antiferromagnetic" orienta- 
tional correlati ons (jl9| ) . Crumpling takes place only if 
1 / Ryfp < 1/ ^/27rr p h and the corresponding persistence 
length is now 



With a computation analogous to the one leading to (|l6| ) 
we obtain the following behaviour: 



r ( r2 
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-arctan(J/i?)ac? 2 , d 2 < - 



87rr p hi?I 

^±-[ln(ad 2 /A)+c(R,I)], ± « d 2 « 1 
d 2 , 



27rr ph 



d 2 » 



27rr„ 



where a = (R 2 + 1 2 ) p and c(R, I) = C + 
[(I 2 - R 2 ) /RI] arctan(//i?) with C = Euler's constant. 

These results show that the model has three possible 
phases. The first is realized when there are no solutions 
to the saddle-point equations in the allowed range of pa- 
rameters. For this choice of parameters, there exist no 
homogeneous, isotropic surfaces. In this phase the sur- 
faces will form inhomogeneous structures. 

If a solution to the saddle-point equations exists, two 
situations can be realized. For large positive stiffness k 
we have R S> /, the asymptotic region is d S> 1/Ry/p 
and / can be neglected. In this region we have 



](d) oc 



1 



y/R~/Pd 



-Ryff>d 



(18) 
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2R AR 2 



(21) 



which is much larger than l/y/fph- Otherwise, the oscil- 
lating superstructure goes over directly into the tension 
dominated region. In this case d E scales logarithmically 
with d for 1/ lyfp y / 2lrr p "h, and p is large not 

because of crumpling but because of the oscillating su- 
perstructure. Indeed there are strong orientational cor- 
relations in this region. Note that the oscillations repre- 
sent a disordered superstructure caused by fluctuations 
on an otherwise homogeneous and isotropic ground-state 
described by the solution of the saddle-point equations. 

4. One might imagine that our disordered superstruc- 
ture undergoes a transition to a crystalline egg-carton- 
type structure || when R — > 0, so that the spectrum of 
transverse fluctuations develops an instability at a finite 
value p = I. However this is not so, as can be seen from 
the explicit expression for p obtained from (KJ), 



P 



16irRI 



it I 2 -R 2 
_ +arctan __ 



(22) 



When lowering R at fixed I one hits a pole where p di- 
verges. This means that the surface crumples before one 
reaches the crystal instability. 
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For symmetry reasons, the transition from the stiff to 
the disordered superstructure phase occurs on the line 
R = I, where k — and the the kernel K develops its 
minimum. 

The fourth-order bending elasticity term dominates 
the fluctuations of surfaces which have both vanishing 
bare tension and stiffness. These can be studied further 
analytically since, for k = (R = I), the saddle-point 
equations become polynomial, with solution 



A = 



128 



1 + 256- 



This gives 



'•ph 



G4 



1 



1 

2ci 



(23) 



(24) 



where a is the following function of the dimensionless 
parameter r/m: 



a 2 = 



1 + 128 r/m + y/l + 256 r/m 



This is the previously announced result. The fourth- 
order bending elasticity term, although irrelevant in per- 
turbation theory, becomes relevant non-perturbatively by 
contributing a term m/64 to the physical surface tension. 
This is the reason why, contrary to stiff membranes j| , 
floppy membranes do not crumple. The physical tension 
can be decreased arbirarily with p remaining in the range 
1 < P < 2. In other words one can lower the two scales 
r and m so that the range of orientational correlations 
d = 1/Ry/p = (4/ayfrn)y/a — 1/2 is always of the same 
order or larger than the inverse of the square root of the 
physical surface tension 1 / ^r p h = 8/ ayfm. 

The corresponding 7-functions are easily obtained as 



/x . d , r 

, . d m 

7 (m) S -A - In - = 2 



(26) 
(27) 



showing the absence of anomalous dimensions for k = 0. 
Correspondingly, we have 



\m 



' A JK £ = £[7(r)-7(m)]=a, (28) 
aA m m 



which means that r/m, and thus also a, are renormaliza- 
tion group invariants. 

5. In conclusion we see that the physics of floppy 
membranes is governed by an infrared-stable fixed-point 
characterized by vanishing tension and a dimension- 
less renormalization group invariant parameter a* = 
lim r _>o,m— ►() o(r, m). At this point, the surface exhibits 
long-range order in which the diagonal elements of the 
correlation functions g a b(£ — £') do not depend on the 
distance, g(d) = Air 2 /a*, and in which the length ( |l7j ) 
scales with the distance in coordinate space like d% — 



(7r 2 p*/a*) d 2 , from which we deduce a Hausdorff dimen- 
sion Dh = 2 for floppy membranes. 
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